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ft ^ T »p»r «ihnws ha t a simple two-link planar manipulator exhibits a phenomenon 
of global Instability in a subspace of its configuration space. A numerical example, as well 
as results of a graphic simulation, is given. 

I. Introduction. 

The problem of unpredictability in deterministic mechanical systems without random 
disturbances whs posed more than a century ago in connection with turbulence motion. A new 
interest in the problem was aroused only recently when it became clear that even low order 
deterministic dynal*lcal systems can be unpredictable from any practical point of view. A 
classical example of aNt^non-Hei les system (1) 12] 
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represents a case of a well posed deterministic conservative system with only two degrees of 
freedom. While (1) cannot be solvetKanalytically for arbitrary initial values, it can be 
integrated numerically. Henon and Heil^ did it in 1964 and the results showed that for the 
system energy above E»l/6, a phase portrait looks seemingly random, while for the E 1/12 <for 
energy levels above E*l/8, while for E<l/8 system (1) exhibits traditional smooth curves 
(Fig. 1). During the last twenty years, existence of chaos in nonlinear dynamic systems 

became a well established fact. 

Another type of chaos that is important for us chn^be well demonstrated by equation <2). 

X{ n+1 ) - 2*X(n) MOD(l) {2) 

For a binary representation of X(n), it simply states that on each step we are shifting the 
decimal point in X one bit right and throwing out an integer part. It is easy to see that (2) 
has an analytical solution (3) 

X(n) - 2 n *X ( 0 ) MOD ( 1 ) <3) 

but to compute a result for a given N. one must know exactly N bits in the initial value of 
X ( 0 ) . 

This example demonstrates a case of orbital instability that was first studied in [1]. 
Equation (3) is a case in which an initial separation between two close solutions grows 
exponentially along the trajectories. Of course for any real dynamical system, orbital 
instability can exist only for the general coordinates that don't increase ths system s total 
energy . It is worth mentioning here that while an exponential "explosion of solutions ^is 
very well known in the Linear System Theory, there it only means that a linear system 
description cannot be used, and the system moves toward its limited circle or breaks down as 
a result of too high stresses. Conversely, an orbital instability in nonlinear systems does 
not lead to an alternative stable equilibrium, and the system description is done, in the 
framework of Newtonian mechanics, without any simplification and linearization. 

IZ. Geometric Approach. 

For our future discussion, a geometrical representation of orbital instability may be 
useful. Let us start from an example of an inertial motion of a single particle M on a 
smooth surface S. In the absence of external forces, a point mass M would me ve along the 
geodesics line on this surface. It is shown in differential geometry [5] that the distance 
between two initially close geodesics is 
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where d„ is an Initial separation, G - Gaussian curvature, and t is a trajectory parameter 
not necessarily time. One can see from (4) that for a negative surface curvature the 
separation increases exponentially. Such a case is shown in Fig. 2. Alternatively, for th 
surfaces with a positive curvature, separation is bounded by its initial value. 


This geometrical representation is very important for a 
stability. Indeed if it is possible to find a space where 
described as an inertial motion of a single particle, one need 
space Gaussian curvature without solving the equations themselves 


question of systea orbital 
the system behavior can be 
examine only the sign qf a 
(*]• 


For the case of conservative f lnlte-degree-of-freedom systems, such a space is vary well 
known. It is a configuration space with a metric tensor corresponding to the structure of the 
system kinetic energy [4]. Let q 1 (i=l..N] denote generalized coordinates, and the kinetic 


system 
energy is 


E - a^q 1 ^ 




then a , , should be used as a new metric tensor 
thereafter, summation is assumed upon repeated 
solution for the free motion of the original system 
a single particle of a unit mass along the geodesic 


for the configuration space. In (5) and 
indexes. In such constructed space, the 
will correspond to an inertial motion of 
lines [5,6]. 


For this metrics triangle, equality is not true any more. Now an elementary arc is 


ds 2 * a^jdq^dq^ 

2 

and ds 2 can be less, greater, or equal to the sum of dq A . 
curvature is connected to this relation. An illustration 
in Fig. 3. 


( 6 ) 

The sign of the resulting Gaussian 
of a two dimensional case is shown 


In the rest of the paper, we are going to show that a free frictionless 
simple mechanical system of a two-link planar manipulator, in the absence of 
demonstrate orbital instability that can be characterized as a weak chaos [4]. 


motion of a 
gravity, can 


III. Solution for a Two-link Arm. 


A model for a two-link planar manipulator is shown in Fig. 4 
be chosen as generalized coordinates q x . q 2 • System kinetic enerqy is 
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Angles fj and f 2 can 
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a 12 * mM*lg*cos( f 2 -f ^ ) (7) 

a 2 2 3 l 2 

where I, and I 2 are moments of inertia, m - mass of the link 2, L - length of link 1, and lg 
is the distance from B to the center of inertia of the link 2. 


The curvature of the resulting two dimensional space can be computed explicitly 

l*a 2 = m*L*lg*(a 2 -f[m*L*lg*sin(f 2 -f 1 )] 2 ) 2 *cos( f 2 -f t ) (8> 


and differential equations for fj and f 2 can be solved numerically. 


One can see from (3) that the sign of the curvature G depends 
Fig. 5 gives the region on (fj,f 2 ) plane where G is positive and 
orbitally unstable. In other words folded arm configurations are 
extended arm configurations are orbitally stable. 


only on the Cbs('f 2 -f 1 ). 
therefore the system is 
orbitally unstable, and 


IV. Numerical Simulation. 

For a numerical simulation, we chose parameter values that had been used for a rea 
manipulator [ 7 ] . 


I. = 0.126, I 2 = 0.075, m = 4.978 
L * 0.27, ig =* 0.0485 


The system was exercised in the following way. For various initial conditions, two arr. z 
were run with a slight difference in their start points. The time history, connquration 
plane trajectories, as well as graphical animation of the arms themselves had been displayed. 
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Alsu we computed and graphically displayed a running estimation of a Lyapunov exponent for 
the f 2 ■ f i 

L(t) * [Ln(d(t)/d 0 )l/t (10> 

where d(t) is the difference between the arms in the value of an internal angle ff 2 (t> — 
fj(t)). It had been shown (8] that a motion is chaotic If 

L(t) — > c >0; when t — > inf (ll) 

Integration was done by using Runge-Kutta with adaptive size steps. To avoid the influence of 
numerical errors, integration was repeated with a tolerance parameter varying more than an 
order of magnitude. All runs gave the salt** • '.ults within desired tolerance. 


V. Results. 

While the fact of positive but not constant curvature over the whole space is a 
necessary but not sufficient condition for the orbital stability, the opposite is true. ■ A 
system is orbitally unstable if the curvature of the space defined by {5), (6) is negative in 
all points. 


Since there are "good" and "bad" regions on the f^ vs f 2 plane {as shown on Fig. 5), to 
get definite results it would be desirable to find an initial condition that would keep the 
system only in the region with positive or negative curvature. However it is clear, that if 
one can hope to find a trajectory that stays completely in the "good" region (G>0), there is 
no trajectory that would stay in the region with orbital instability. Indeed, any solution 
f 2 - fi Is of M0D(2PI), and the separation can not grow indefinitely without forcing an 
"unfold" . 


for the 
arm to 


In our simulation we found a case when an arm starting from unfolded position would stay 
there. In that case the original difference between the arm almost did not grow as can be 
seen from Fig. 6a. This case will be further referred to as Case I. Its initial conditions 
were : 


arml : f, = 80, f, = 95 (dg], = 300. f 2 - 0 [dg/sec] 

arm2 : fj = 80. f 2 = 97 (dg]. t l = 300. f 2 = 0 (dg/sec] 

For the arm started from a folded position a small initial difference grew very fast as 
could be seen from Fig. 6b. Fig. 6c shows arms in one of the intermediate positions. It is 
how far apart they become. The initial conditions for the Case II were 

arml : f, * 80. f 2 = 200 (dgj, = 300. f 2 3 0 (dg/sec] 

arm2 : fj = 80. f 2 ■= 202 ( dg ] . f i - 300. f 2 = 0 [dg/sec] 

The estimate of the Lyapunov exponent for the Case I clearly goes to zero, while 
Case II stayed positive at least during the time of observation - Fig. 7a, b. It is not clear 
nevertheless that it will never go to zero, due to the reasons described above (crossing both 
good as well as bad regions) . 


VI. Discussion. 

A phenomenon of chaotic motion has been theoretical ly found and numerical ly illustrated 
for a simple mechanical system of a two-link manipulator. It nas been shown that a * ** e 
arm is oroitally unstable, opposite to an extended one. 

While the assumption of a nonfriction, zero gravity environment is quite unrealistic, we 
believe that our finding warrants further and more detailed investigations of the described 
phenomena. It is quite possible that low-friction, many degrees of freedom redundant flexible 
arms of the future will exhibit more complicated behavior that could lead to orbital 
instability under more realistic conditions. 

The importance of the geodesics in the configuration space with metric a A j as minimum 
time trajectories has been recently shown by K.G. Shin and N.D. McKay (9], Our result 
suggests that an open loop control should not be used in the region with negative curvature 
since trajectories there rapidly diverge. 

A question of a closed loop control in the area of negative curvature also deserves 
more detailed investigation. While it had been shown that a simple PID control under the same 
conditions (no friction and no gravity) makes an arbitrary robot manipulator asymptotically 
stable [101 it is not clear how trajectory curvature affects an admissible sampling rate. 
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VII. Conclusion. 

Using a geometric approach we have shown that a simple robot manipulator can be orbi tally 
unstable depending on Its configuration. A numerical simulation supported this finding. Me 
feel that further efforts in this direction will help better understanding of the dynamical 
properties of such complicated nonlinear systems as robot manipulators. 


Acknowl edgement 

This work was done at the Jet Propulsion Laboratory, California Institute of Technology under 
a contract with the National Aeronautics and Space Administration (NASA). 


Literature : 

Cl] W. Thompson, P.G .Tout , "Treatise of Natural Philosophy" voli, 
part 1, p.416 1879. 

[2] M.Henon, C.Heiles, Astron.J. 69,73 . 196* 

[3] Joseph Ford "How random is a coin toss?" 

Physics Today, April 1983, pp 40-47. 

[4] M. Zak "Two Types of Chaos in Non-Linear Mechanics" 

X. J of Nonlinear Mechanics, v. 20, Nov 4, pp. 297-308, 1985. 

[5] J.L. Singe "On the Geometry of Dynamics". Phil. Trans. R. 

Soc. Lund., Ser A, 226, pp . 31-106 (1096) 

[6] A. I. Lurie "Analytical Mechanics" (Russian) Moscow, 1961 

[7] 0. Sato.H. Shimojlma, Y.Kltamura "Minimal Time Control of a 
manipulator with 2 degrees of freedom". Bulletin JSME vol 
26, no. 218, Aug. 1983 

[8] A . J . Lichtenberg, M.A. Lieberman 

"Regular and Stochastic Motion", Springer-Verlag, NY, 

1983 

[9] K.G.Shin, N.D. McKay "Selection of Near-minimum Time 

Geometric Paths for Robotic Manipulators" 

IEEE Tr. AC v. AC-31 No. 6, June 1986. pp. 501-511 

[10] M.Takegaki, S.Ariraoto "A new Feedback Method for Dynamic 
Control of Robot Manipulators". JDSMC June 1981, v.102, pp. 119-125 



v 



figure 1 


FIGURE 2 


232 



THETOZ s/s TIME ' IHETftZ s/% TIMI 



